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BACKGROUND 

20 

The ability to measure the three-dimensional structure of 
biological tissues is important, but common methods used to 
leave out basic information about the molecular composition 
and metabolic behavior of the tissue imaged. This molecular 
composition and metabolic behavior information may yield 25 
valuable scientific data on the behavior of biological systems, 
and would be of great clinical diagnostic value for finding 
diseases such as cancer. Much of the focus of biological and 
medical imaging today is to gain information about compo- 
sition. 30 

Functional magnetic resonance imaging (FMRI) utilizes 
contrast agents specific to particular molecular species or 
metabolic processes to provide specific information about the 
three-dimensional (3-D) location of these species or pro- 
cesses. This method is quite versatile, and poses relatively 35 
little risk to the patient, but is limited in practice to a resolu- 
tion over 100 microns. Positron Emission Tomography (PET) 
and Single Photon Emission Computed Tomography 
(SPECT) utilizes a radiolabeled metabolic molecule to pro- 
vide 3-D measurements of utilization of these molecules by 40 
sensing their gamma -ray emissions. Unfortunately, these 
methods have been cost prohibitive as diagnostic tools, and 
expose the patient to ionizing radiation. Fluorescence micros- 
copy labels relevant biological structures and processes with 
a fluorescent compound which can be measured by conven- 45 
tional microscopy or confocal microscopy. The two-photon 
variant uses ultrafast pulses so that emission and excitation 
frequencies can be clearly separated, and also utilizes the 
additional selectivity that the intensity-sensitive nature of 
two-photon excitation provides. While some biological struc- 50 
tures produce natural fluorescence, in general externally 
introduced fluorescence markers must be used, which can 
interfere with biological processes and are often toxic. 

Molecules frequently have molecular resonance frequen- 
cies that are due to the electromagnetic attractions of atoms in 55 
the molecule. These frequencies are those of molecular vibra- 
tions, molecular rotational motions, the excitation of elec- 
trons to higher energy states, and occasionally finer structures 
such as hyperfine interactions and optical-magnetic proper- 
ties. These properties are present without the introduction of 60 
any external contrast molecule. These frequencies are usually 
in the mid-infrared, corresponding to photons of 1.5-50 
microns of wavelength. Unfortunately, they cannot be 
directly excited by electromagnetic radiation of the same 
frequency because when they are in tissue, the surrounding 65 
water absorbs almost all of these frequencies. The range of 
wavelengths that the tissue is relatively transparent to is 0.6- 


2 

1 .5 microns. Therefore multiphoton nonlinear processes need 
to be employed to probe these resonances. The photons to 
stimulate and record the processes are typically in a region 
where the tissue is not absorbing, so that they can reach the 
tissue feature and be measured from the feature. 

Raman spectroscopy, first discovered in 1928, uses 
molecular resonance features of frequency Aoo to split a pho- 
ton of frequency 00 into another photon of frequency 00 -Aoo 
and a resonance excitation of frequency Aoo. The presence of 
photons at frequency 00 -A 00 identifies the concentration of the 
resonance feature. This process is in practice very weak and 
requires large amounts of power to produce any detectable 
amounts of photons. This weakness is due to the fact that the 
probability of a Raman excitation process to occur is propor- 
tional to the number of photons at frequency 00 -Aoo already 
present, of which there are typically few or none. Since pho- 
tons that would be emitted by Raman excitation at frequency 
00 -A 00 are indistinguishable from the incoming radiation that 
stimulates them, this is not a viable technique for achieving 
molecular sensitivity. 

Coherent Anti- Stokes Raman Scattering (CARS) is 
another nonlinear spectroscopy technique that unlike conven- 
tional Raman spectroscopy, allows all of the photons neces- 
sary to stimulate the process to be introduced into the tissue 
by the illuminating source. This enables the probability of a 
CARS interaction to be increased to a (theoretically arbi- 
trarily) high level so that a sufficient number of photons can 
be produced as to enable detection within a reasonable time 
period. It is essentially two stimulated Raman scattering pro- 
cesses in parallel. Two photons, a “pump” of frequency co 1 
and a “Stokes” of frequency a > 2 illuminate the tissue. They 
must be separated in frequency by co 1 -co 2 =Aco, which is the 
frequency of the molecular resonance. When molecules of the 
target molecular species are present, the resonance will be 
excited, and the pump photon will be converted to the same 
frequency as the Stokes photon. This is the first stimulated 
Raman scattering process. Another photon may arrive at fre- 
quency co 3 that will stimulate the emission of the excitation 
from the resonance, so that the energy of the photon of fre- 
quency oo 4 and the excitation are converted to a new photon of 
frequency o) 4 =(% +A 00 , called the “anti-Stokes” photon. The 
presence of this photon of frequency cg 4 indicates that a 
CARS process has taken place and indeed a molecule with the 
resonance feature is present. Often the “pump” beam is used 
as the photons of frequency o) 3 , so that 0 ) 3 = 0 )! and 0 ) 4 = 20 ^- 
o) 2 . Since the photon of o ) 4 is not the same frequency as one of 
the illuminating photons, and is typically within the transpar- 
ency range of the tissue, it is easily discriminated from the 
incoming radiation. FIG. 1A shows an energy-level diagram 
for CARS, and FIG. IB shows an energy-level diagram for 
Coherent Stokes Raman Scattering. 

CARS microscopy uses the CARS process to look for the 
presence of a molecular species, but does not require any 
foreign substances to be introduced into the tissue. It scans the 
illumination point-by-point through the tissue and measures 
the number of generated anti-Stokes photons. When a three- 
dimensional mesh of points has been scanned, a complete 
three-dimensional picture of molecules of that resonance can 
be shown. Since CARS is a nonlinear process (and therefore 
is intensity sensitive), efficient conversion only occurs at the 
focus of the illumination, which can be made very tight (typi- 
cally a half micron in both the axial and lateral directions). 
Therefore the resolution can be made many orders of magni- 
tude better than MRI, which is the probably the largest com- 
petition for clinical use for similar purposes. Unfortunately, 
the penetration is usually rather low (limited to about 500 
microns). A further shortcoming is that CARS microscopy 
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measures the total number of anti-Stokes photons, or power, 
from the sample. However, the optical field contains temporal 
structure in the phase that is averaged out by power detection 
because photodetector response time is orders of magnitude 
slower than the oscillations of the optical field. The time scale 
on which the optical pulse varies (which is typically picosec- 
onds or femtosecond time scales) is far too fast for photon 
detection equipment or electronics to detect (the fastest of 
which may detect 25 ps time scales). 

Optical coherence tomography (OCT) is an emerging 
high-resolution medical and biological imaging technology 
[15-21]. OCT is analogous to ultrasound B-mode imaging 
except reflections of low-coherence light are detected rather 
than sound. OCT detects changes in the backscattered ampli- 
tude and phase of light. 

Cross-sectional OCT imaging is performed by measuring 
the backscattered intensity of light from structures in tissue. 
This imaging technique is attractive for medical imaging 
because it permits the imaging of tissue microstructure in situ, 
yielding micron-scale imaging resolution without the need 
for excision and histological processing. Because OCT per- 
forms imaging using light, it has a one- to two-order-of- 
magnitude higher spatial resolution than ultrasound and does 
not require contact with tissue. 

OCT was originally developed and demonstrated in oph- 
thalmology for high-resolution tomographic imaging of the 
retina and anterior eye [22-24] . Because the eye is transparent 
and is easily optically accessible, it is well-suited for diag- 
nostic OCT imaging. OCT is promising for the diagnosis of 
retinal disease because it can provide images of retinal pathol- 
ogy with 10 pm resolution, almost one order-of-magnitude 
higher than previously possible using ultrasound. Clinical 
studies have been performed to assess the application of OCT 
for a number of macular diseases [23,24]. OCT is especially 
promising for the diagnosis and monitoring of glaucoma and 
macular edema associated with diabetic retinopathy because 
it permits the quantitative measurement of changes in the 
retinal or retinal nerve fiber layer thickness. Because morpho- 
logical changes often occur before the onset of physical 
symptoms, OCT can provide a powerful approach for the 
early detection of these diseases. 

Recently, OCT has been applied for imaging a wide range 
of nontransparent tissues [16,17,25-27]. In tissues other than 
the eye, the imaging depth is limited by optical attenuation 
due to scattering and absorption. A “biological window” 
exists in tissue where absorption of near-infrared wave- 
lengths is at a minimum and light can penetrate deep into 
highly-scattering tissue (FIG. 15 ) [28]. Because optical scat- 
tering decreases with increasing wavelength, OCT in non- 
transparent tissues has routinely used 1.3 \im wavelength 
light for imaging. In most tissues, imaging depths of 2-3 mm 
can be achieved using a system detection sensitivity of 1 1 0 dB 
(1 part in 10 11 ). OCT has been applied to image arterial 
pathology in vitro and has been shown to differentiate plaque 
morphology with superior resolution to ultrasound [17,29]. 

Imaging studies have also been performed to investigate 
applications in gastroenterology, urology, and neurosurgery 
[30-32]. High resolution OCT using short coherence length, 
short-pulse light sources, has also been demonstrated and 
axial resolutions of less than 5 pm have been achieved [33, 
34]. High-speed OCT at image acquisition rates of 4 to 8 
frames per second for 500 to 250 square pixel images has 
been achieved [35] . OCT has been extended to perform Dop- 
pler imaging of blood flow and birefringence imaging to 
investigate laser intervention [36-38]. Different imaging 
delivery systems including transverse imaging catheters and 
endoscopes, and forward imaging devices have been devel- 
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oped to enable internal body OCT imaging [39,40]. Most 
recently, OCT has been combined with catheter-endoscope- 
based delivery to perform in vivo imaging in animal models 
and human patients [41 -44] . 

5 Apart from medical applications, OCT has been demon- 
strated as an emerging investigational tool for cell and devel- 
opmental biology. OCT has imaged the development of 
numerous animal models including Rana pipiens and Xeno- 
pus laevis (Leopard and African frog), and Brachydanio rerio 
to (zebrafish) [45-46]. High-speed OCT imaging has permitted 
the morphological and functional imaging of the developing 
Xenopus cardiovascular system, including changes in heart 
function following pharmacological interventions [47] . 
High-resolution imaging has permitted the real-time tracking 
1 5 of cell dynamics in living specimens including mesenchymal 
cell mitosis and neural crest cell migration [48]. OCT is 
advantageous in microscopy applications because repeated 
non-invasive imaging of the morphological and functional 
changes in genetically modified animals can be performed 
20 overtime without having to histologically process multiple 
specimens. The high-resolution, cellular- imaging capabili- 
ties suggest that OCT can be used to diagnose and monitor 
early neoplastic changes in humans. 

The ability of OCT to perform optical biopsies, the in situ 
25 imaging of tissue micro structure at near-histological resolu- 
tion, has been used to image morphological differences 
between normal and neoplastic tissue. OCT images of in vitro 
neoplasms of the female reproductive tract [49], the gas- 
trointestinal tract [50], and the brain [51] have been investi- 
30 gated. Optical differences between normal and neoplastic 
tissue were evident, but primarily for late-stage changes. Still, 
situations exists were no inherent optical contrast exists 
between normal and pathologic tissue, such as in early- stage, 
pre-malignant tumors or in tumors which remain optically 
35 similar to normal tissue. 

In the past, OCT has found numerous medical and biologi- 
cal applications. However, the imaging technique has relied 
largely on the inherent optical properties of the tissue to 
provide contrast and differentiate normal from pathological 
40 tissue. Phospholipid-coated perfluorobutane microbubbles 
(ImaRx Pharmaceutical, Tucson, Ariz.) have been used as a 
contrast agent for OCT; although they produce a strong OCT 
signal, blood and tissue also produce a fairly strong OCT 
signal, and the effects of this contrast agent in vivo on the 
45 visualization of blood vessels are subtle. 

SUMMARY OF THE INVENTION 

In a first embodiment, the invention provides a method of 
50 examining a sample, which includes: exposing a reference to 
a first set of electromagnetic radiation, to form a second set of 
electromagnetic radiation scattered from the reference; 
exposing a sample to a third set of electromagnetic radiation 
to form a fourth set of electromagnetic radiation scattered 
55 from the sample; and interfering the second set of electro- 
magnetic radiation and the fourth set of electromagnetic 
radiation. In this embodiment, the first set and the third set of 
electromagnetic radiation are generated from a source. More- 
over, at least a portion of the second set of electromagnetic 
60 radiation is of a frequency different from that of the first set of 
electromagnetic radiation, and at least a portion of the fourth 
set of electromagnetic radiation is of a frequency different 
from that of the third set of electromagnetic radiation. 

In a second embodiment, the invention provides a method 
65 of forming an image of a sample, which includes: exposing a 
reference to a first set of electromagnetic radiation, to form a 
second set of electromagnetic radiation scattered from the 
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reference; exposing a sample to a third set of electromagnetic 
radiation to form a fourth set of electromagnetic radiation 
scattered from the sample; forming a digital data set corre- 
sponding to the sample; and converting the data set into an 
image. In this embodiment, the data set is formed by interfer- 
ing the second set of electromagnetic radiation and the fourth 
set of electromagnetic radiation. Also, the first set and the 
third set of electromagnetic radiation are generated from a 
source. Moreover, at least a portion of the second set of 
electromagnetic radiation is of a frequency different from that 
of the first set of electromagnetic radiation, and at least a 
portion of the fourth set of electromagnetic radiation is of a 
frequency different from that of the third set of electromag- 
netic radiation. 

In a third embodiment, the invention provides a device for 
examining of a sample having an oscillator, a reference gen- 
erator that is optically coupled to the oscillator, a microscope 
that is optically coupled to the oscillator, a demodulator that 
is optically coupled to the reference generator and the micro- 
scope, and a recorder that is coupled to the demodulator. 

In a fourth embodiment, the invention provides a method of 
examining a sample that includes exposing a sample to a first 
set of electromagnetic radiation to form a second set of elec- 
tromagnetic radiation scattered from the sample, and inter- 
fering the second set of electromagnetic radiation with a third 
set of electromagnetic radiation. In this embodiment, the third 
set of electromagnetic radiation is phase-coherent with the 
first set of electromagnetic radiation, at least a first portion of 
the second set of electromagnetic radiation is of a frequency 
different from that of the first set of electromagnetic radiation, 
and at least a portion of the third set of electromagnetic 
radiation is of the same frequency as the first portion of the 
second set of electromagnetic radiation. 

In a fifth embodiment, the invention provides a method of 
forming an image of a sample that includes exposing a sample 
to a first set of electromagnetic radiation to form a second set 
of electromagnetic radiation scattered from the sample, form- 
ing a digital data set corresponding to the sample, and con- 
verting the data set into an image. In this embodiment, the 
forming of the image includes interfering the second set of 
electromagnetic radiation and a third set of electromagnetic 
radiation. Moreover, the third set of electromagnetic radiation 
is phase-coherent with the first set of electromagnetic radia- 
tion. In addition, at least a first portion of the second set of 
electromagnetic radiation is of a frequency different from that 
of the first set of electromagnetic radiation, and at least a 
portion of the third set of electromagnetic radiation is of the 
same frequency as the first portion of the second set of elec- 
tromagnetic radiation. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1. Coherent Anti-Stokes Raman Scattering and 
Coherent Stokes Raman Scattering energy -level diagrams. 

FIG. 2. Basic block diagram of NIVI 

FIG. 3. Example laser configurations are shown that pro- 
duce two pulses of frequencies and oo 2 overlapped. We 
assume that the delays have been set correctly to overlap 
them. 

FIG. 4. Methods of shaping a broadband pulse into a pulse 
with beat frequency Aoo. 

FIG. 5. Reference Signal Generator Implementation 

FIG. 6. Configurations for full field CARS. 

FIG. 7. Translated serial-point scanning configurations 

FIG. 8. Beam-steered serial-port scanning configurations 

FIG. 9. Temporal-ranging based NIVI 

FIG. 10. Full field cross-correlator demodulator 
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FIG. 11. Temporal cross -correlator for a serial -point scan- 
ning microscope. 

FIG. 12. Shown are two configurations that utilize linear 
photodetector arrays to measure multiple samples of the 
5 cross -correlation simultaneously. 

FIG. 13. Diagram of an implementation of the recorder. 

FIG. 14. Electrical cross-correlation signal quadrature 
amplitude demodulator. 

FIG. 15. “Biological window” in tissue where absorption 
to of near-infrared wavelengths is at a minimum and light can 
penetrate deep into highly-scattering tissue. 

DETAILED DESCRIPTION 

15 Nonlinear interferometric vibrational imaging (NIVI) is a 
method used to measure the three-dimensional distribution of 
molecular species in various samples (biological or other- 
wise). Its basic operation is to stimulate the excitation of 
molecular bonds with particular resonance frequencies, and 
20 then use these excitations to produce radiation distinct from 
the excitation that can be measured. The physical process of 
excitation and stimulation of radiation is called Coherent 
Anti-Stokes Raman Scattering (CARS). Unlike previous 
methods that use CARS in microscopy to probe for the pres- 
25 ence of molecular species, NIVI utilizes a heterodyne 
approach where a reference signal is separately generated and 
interferomerically compared to the signal received from the 
sample, allowing the signal to be formed into an image in the 
same way as OCT. In this way, additional information can be 
30 inferred from the emitted radiation such as the distance to the 
sample and phase information that yields additional structure 
of the molecular bonds. It also has other advantages in sen- 
sitivity and the ability to screen out background radiation that 
is not produced by the sample. It also can allow more flex- 
35 ibility in the choice of laser illumination source, because the 
coherent detection process does not rely on photon frequency 
alone to discriminate emitted radiation. 

The term “scattered photon” means photons scattered by a 
sample, which include linearly scattered photons and non- 
40 linearly scattered photons, such as CARS and CSRS photons . 

The term “image” means data produced by receipt of elec- 
tromagnetic radiation, which may or may not be formed into 
a picture viewable by the human eye. This includes images 
produced directly onto a medium such as film or video. 

45 The phrase “frequency range of infra-red to ultraviolet” 
means electromagnetic radiation having a frequency of 10 12 
to 10 17 Hz, which excludes radio waves, microwaves, X-rays 
and gamma rays. The term “light” means visible light. 

The method of NIVI was developed at first because CARS 
50 microscopy based on non-interferometric detection is incom- 
patible with the scanning modes of OCT. When the method of 
CARS microscopy is employed, information about the echo 
time of the optical signal is lost because the photodetector is 
not nearly fast enough to respond to the relatively instanta- 
55 neous return of the return pulse. In OCT a reference wave is 
used so that the relative time delay between the reference 
wave and the returned echo wave can be determined. Nonlin- 
ear methods (such as conventional Raman spectroscopy) are 
unsuited to integration with OCT because the emitted signal 
60 is incoherent with respect to the excitation. However, CARS 
processes preserve coherence, and therefore enable interfero- 
metric methods to be used. By using a nonlinearly generated 
reference wave, the time of arrival of a CARS signal returning 
from the tissue can be interferometrically compared to the 
65 arrival time of the reference signal . In this way, the methods of 
OCT and CARS imaging can be integrated. In addition, in a 
manner analogous to how OCT can measure the dispersion of 
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scatterers in the medium, NIVI can measure the dispersion of 
the response of molecules to the excitation radiation. This is 
because NIVI measures the relative phase between the refer- 
ence signal and the sample signal. This dispersion should 
contain information about the resonance structure of the mol- 5 
ecule over and above what can be measured using CARS 
microscopy. 

The Nonlinear Interferometric Vibrational Imaging 

Method io 

The purpose of NIVI is to measure the temporal field 
radiated by samples that are stimulated using the CARS tech- 
nique. When the field is scanned through the sample, an 
image can be formed of the molecular contents of the sample. 15 
As illustrated in FIG. 2, a preferred embodiment NIVI 200 
has the following components: 

Oscillator. This oscillator 201 produces an optical field that 
can excite the resonance mode of the target molecule through 
a nonlinear technique (usually stimulated Raman scattering), 20 
and also the photon to stimulate the photon that is measured 
(also usually through stimulated Raman scattering). The 
combination of these processes is called CARS. 

Reference Generator. The reference signal generator 203, 
which can sometimes be incorporated into the oscillator, con- 25 
verts part of the oscillator signal to a reference signal that can 
be used in the demodulator section. It acts as a known signal 
that demodulates the unknown signal from the sample in the 
interferometer. 

Microscope. The microscope 205 delivers the field pro- 30 
duced by the oscillator to the sample, and collects the field 
emitted by the sample. The excitation field is usually deliv- 
ered by a microscope objective, where the oscillator field is 
focused tightly or sparsely, depending on the scanning 
method. This focus is scanned through the tissue, and based 35 
on the signal received from each tissue volume an image can 
be formed. When the oscillator signal is delivered to the tissue 
containing a molecule with a compatible resonance, a non- 
linear process such as CARS can occur and produce a new 
sample signal (called the “anti-Stokes” for CARS processes). 40 
This sample signal serves as an indicator of the presence of 
the molecular resonance, and also provides additional infor- 
mation about the molecule through the temporal structure of 
the sample signal. 

Demodulator. The demodulator 207 combines the signal 45 
received from the sample with the reference signal. This is 
typically achieved by constructing an interferometric cross- 
correlator. The cross-correlation of the two signals is then 
measured by a single photodetector or array of photodetec- 
tors. The power received by these photodetectors allows the 50 
cross -correlation signal to be inferred, from which the tem- 
poral response signal from the sample can be also inferred. 
With knowledge of the physics of the molecule, the presence 
of and potentially properties of the molecule being tested can 
be inferred from its temporal response. 55 

Recorder. The data recorder 209 records the data measured 
by the demodulator. This data can be digitally processed to 
produce an image that a human operator can interpret. 

Each of these modules can be implemented in a variety of 
different ways that can be tailored to various data acquisition 60 
needs. In addition, while these units are the basic units of the 
invention, often the parts can be consolidated to simplify 
implementation or reduce cost. While the basic block design 
could be construed as that of a standard interferometric 
microscope, nonlinear processes are occurring in the refer- 65 
ence generator and “Microscope” sections that allow the reso- 
nance information of the sample to be pumped. 
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Each of these units will be detailed presently. 

1 . Oscillator. 

The oscillator produces the electromagnetic field that 
stimulates the excitation of the resonance to be probed. It also 
provides the photon that stimulates the output photon that is 
detected as evidence of CARS or CSRS. There are many 
types of oscillators and fields that can produce CARS. Each 
pulse produced by the oscillator should be nearly identical so 
that it can excite consistent signals in the reference generator 
and sample. If the oscillator produces too variable of a signal, 
the signals from the reference generator and sample may 
change and produce signals that can be confused with noise 
sources. Variability in the oscillator output is a noise source in 
itself that adds uncertainty to what the expected demodulated 
signal should be. 

The conventional way to produce CARS is to send in two 
overlapped optical pulses, one of which at frequency co l5 the 
pump, and the other at co 2 , the Stokes pulse, where 
o)i— a) 2 = Aa), where Aco is the resonance frequency of the 
molecules of interest. These pulses produce a beat frequency 
of Aco that manifests itself in the magnitude of the optical 
field. In linear time-independent optics, systems are sensitive 
only to the frequencies of the optical pulses themselves, and 
not any beats they may produce together. However, with 
sufficient intensity the intensity envelope may itself stimulate 
the molecule. By choosing two pulses that produce beats of 
this frequency, we can stimulate the molecule with two wave- 
lengths that the tissue is transparent to. Once the resonance is 
stimulated, another photon of frequency oq (in CARS), or of 
frequency oo 2 (in CSRS) stimulates the emission of a fourth 
photon, which is of frequency 2<j0i -ecu for CARS, and 2oo 9 - 
uq for CSRS. 

Systems that can be used to produce these two frequencies 
are shown below in FIG. 3. A common configuration to pro- 
duce pulses of these two wavelengths that are overlapped in 
time is to have a pulsed laser produce one of these wave- 
lengths, split off of some of its energy, and use this energy to 
produce a second pulse of a lower or higher frequency. In one 
configuration 300, a pump laser 301 pumps a dye laser 303, 
for example a doubled Nd: YAG pump laser at 532 nm pump- 
ing a tunable dye laser. In another possible configuration 302, 
a pump laser such as a Ti-sapphire oscillator pumps an optical 
parametric oscillator (a device that converts pulses to lower 
frequencies) 305. 

In yet another configuration 304, the pump laser pumps a 
regenerative amplifier 307, such as a Ti-sapphire regenerative 
amplifier. The regenerative amplifier then pumps an optical 
parametric amplifier (another frequency conversion device) 
309. Alternatively, as illustrated in configuration 306, the 
pulses of each wavelength are generated by two separate 
pump lasers, and the time overlap is maintained by a circuit 
311 that synchronizes the two sources. In another configura- 
tion 308, the pump laser pumps a continuum light generator 
3 13 , generating broadband light which is filtered by a filter for 
the two wavelengths with group velocity dispersion correc- 
tion 315. 

While directly generating the two frequencies and super- 
imposing them to produce beats is the most common way to 
stimulate CARS, this method has some disadvantages for the 
method of NIVI. In CARS and CSRS, there are two types of 
generated signals. Resonant signals depend on the presence 
of a molecule of a particular resonance frequency to be 
present to generate the CARS signal. Another component, 
nonresonant CARS, does not require a particular frequency to 
perform conversion. Nonresonant CARS depends on the peak 
intensity in the signal, while the resonant component can 
build up from many beat periods and so therefore can be 
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spread out in time. Because of this, it is advantageous to 
spread the CARS signals in time to reduce the nonresonant 
component. 

However, when the two signals are discretely generated 
and are transform-limited (are not chirped in time), the only 
way to broaden the signals in time is to reduce their band- 
width. To achieve sufficient power- spectral-density to cause 
efficient conversion, the pulses must either generated by a 
low-bandwidth laser, or much power will be wasted in filter- 
ing a higher bandwidth signal. Unfortunately, the range reso- 
lution in OCT, when temporal ranging is used, is inversely 
proportional to the illumination bandwidth. This requirement 
for high bandwidth conflicts with the requirement for small 
bandwidth for resonance specificity. It would be desirable to 
come up with an alternate oscillator configuration that would 
preserve the resonance specificity of the low bandwidth 
pulses, but actually utilize high bandwidth signals. 

Since the nonlinear excitation of the resonant molecule 
depends more on the beats produced than on the bandwidth 
used to produce them, it would be desirable to take a broad- 
band pulse and reshape it into a signal with the required beat 
frequency. Recent advancements have made pulsed sources 
of very large bandwidth. Some of the methods to do this are 
high-bandwidth Ti-sapphire oscillators, dispersion compen- 
sated mirror Ti-sapphire oscillators, double chirped-mirror 
Ti-sapphire oscillators, and continuum generation sources. 
The optical field produced by these sources can be shaped into 
a field with the beats at the required frequency. 

One such method that has been demonstrated in the litera- 
ture is shown in FIG. 4. A source of laser pulses from a laser 
source 401 is sent into a Fourier-plane pulse shaper 403 that 
utilizes a spatial-light-modulator (e.g. liquid crystal modula- 
tor or acousto-optic modulator). The Fourier-plane pulse 
shaper enables each frequency in the pulse to have its phase 
and/or amplitude altered. By applying the correct phase and 
amplitude to each incoming frequency, the incoming signal 
can be convolved with an essentially arbitrary signal. The 
pulse shaper is set up to reshape the incoming pulse by apply- 
ing a period phase or amplitude perturbation in the Fourier 
domain with period Aco/N, where N is a positive integer. This 
will transform a single pulse into a train of pulses that are 
separated in time by 2jt/Aco. If only a phase perturbation is 
used, the power of the signal can be maximally preserved. The 
larger the integer N is, the longer the pulse train will be, and 
none of the bandwidth of the original pulse will be lost. 
However, most pulse shapers have a limited number of con- 
trollable frequencies, limiting the practical size of N. 

One advantage of spatial-light-modulator based pulse 
shapers is that there is typically a wide range of pulse shapes 
that can be achieved, and the spatial-light-modulator can 
often be controlled automatically by a computer. The com- 
puter can then adjust the spatial-light-modulator to achieve 
maximum signal from the sample in a feedback loop. This 
may allow automatic correction of dispersion or aberrations 
introduced by the optics of the system, and will permit more 
flexibility in probing the molecule because of the tunability of 
the pulse shapes. 

The pulse-shaper in the “Pulse-shaper type NIVI oscilla- 
tor” 400 is well described in the literature. It consists of two 
diffraction gratings, which disperse and recombine the fre- 
quencies, two lenses that focus each frequency to a point and 
recollimate each frequency, and a pulse shaper placed at the 
focal plane to separately operate on each frequency. The pulse 
shape is altered by dispersing each frequency to a separate 
angle, and then imaging each frequency to a separate point on 
the spatial -light -modulator. Alternatively, an etalon may be 
used to shape the amplitude of the pulse periodically. Unfor- 
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tunately, while this would be simpler, it modifies the spectrum 
of the pulse and therefore introduces artifacts into the NIVI 
image. 

An alternative method is to take a pulse and impart a linear 
5 chirp to it. A linear chirp turns a pulse into one where the 
frequency rises or falls at a linear rate as a function of time. 
This rate is characterized by a constant a, which is the change 
in frequency per unit time. It is called “chirped” because of 
the noise of the equivalent sound wave. If two copies of the 
to chirped pulse are created, delayed with respect to each other 
by imparting a variable time delay 407, and recombined, the 
resulting pulse will have two simultaneous frequencies that 
will rise or fall together at the same linear rate, but always be 
separated at a given instant by the same frequency. If this 
1 5 separation frequency is chosen to be Aoo, then the envelope of 
the pulse will be modulated by beats of this frequency. This 
method is especially convenient because the probed reso- 
nance frequency can be adjusted easily by adjusting the time 
delay between the two chirped pulses, which is relatively easy 
20 and inexpensive. This will enable a NIVI instrument that can 
be rapidly and easily adjusted to scan a wide range of molecu- 
lar resonances. Systems based on tunable frequency sources 
will likely be much more difficult to dynamically change 
reliably and often. 

25 The ability to linearly chirp a pulse is well known in the 
literature. It can be accomplished with a pulse shaper 405 
having a combination of prisms, diffraction gratings, lenses, 
mirrors, and/or dispersive materials. Combinations may be 
required to ensure that the resulting chirp is linear and does 
30 not contain significant amounts of higher-order dispersion. 
Higher-order dispersion would limit the resolution to which 
the resonance could be addressed and exclude other nearby 
frequency resonances. In a typical Chirped CARS NIVI setup 
402, the chirp rate required would be fixed and the chirp rate 
35 should need little or no adjustment in the field. Measuring 
devices such as Frequency Resolved Optical Gating can test 
whether a chirped pulse is linearly chirped. 

When using high-bandwidth excitation for CARS, it is 
important to filter out the entire bandwidth of excitation 
40 before detection so it does not interfere with detection of the 
emitted anti-Stokes light (Stokes for CSRS), because the 
nonlinear emission can not be easily distinguished from the 
much larger linearly scattered excitation light. However, this 
linearly scattered light contains the same structure that con- 
45 ventional OCT imaging does, and may be used to measure 
this information at the same time that a NIVI image is 
recorded. This will be convenient for superimposing OCT 
and NIVI data onto the same image, because the acquisition 
of both types of data can be designed into the same instru- 
50 ment. This can be implemented in practice by using a dichroic 
beamsplitter to separate the excitation and response radiation, 
and measuring the cross-correlation of the two frequency 
bands separately using a cross -correlation demodulator. 

With high-bandwidth sources where the entire bandwidth 
55 need not be utilized to produce the excitation field, it is 
possible to use the upper end (for anti-Stokes) or lower end 
(for Stokes) of this bandwidth as a reference field, eliminating 
the need to separately generate a reference field. However, the 
frequencies of the reference must occupy the same band as the 
60 received CARS/CSRS light from the sample. Some sources, 
especially continuum generation sources, will likely generate 
much more bandwidth than needed for pulse shaping and 
therefore will probably have this extra bandwidth available 
for this use. While a different process from CARS/CSRS 
65 typically generates this light, it will likely remain phase- 
coherent with the CARS/CSRS light and therefore should be 
useful as a reference. Phase-coherence depends on the 
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mechanism of pulse/continuum generation and therefore its 
phase-coherence stability properties of a particular source 
type must be established before it is suitable for this purpose. 
A dichroic beamsplitter may be used to separate the fre- 
quency band corresponding to the response radiation from the 
oscillator energy, so that it may be utilized as a reference 
signal. 

It is also possible to simultaneously stimulate the excita- 
tion of several resonances if the sample is illuminated with the 
appropriate pump and Stokes/anti- Stokes beams. For 
example, in CARS a narrowband pump signal and a wideband 
Stokes can be used to address many resonances simulta- 
neously. This is called multiplex CARS and can be extended 
to CSRS with a broad anti-Stokes wavelength range. This 
may be used to measure the presence of several molecular 
resonances simultaneously in the sample. In addition, if sev- 
eral excitations can be produced in the same molecule simul- 
taneously, the molecule will evolve to various quantum states 
depending on the relative amplitude and timing between the 
CARS/CSRS stimulating signals for each resonance. This 
may be produced by sending in multiple pairs of Stokes/anti- 
Stokes wavelengths and pump beams in with varying time 
delays between them. By varying the time delay between 
excitations, the molecule can be made to prefer Stokes or 
anti-Stokes emissions from a particular quantum state. This 
way, the amount of anti-Stokes radiation generated from each 
quantum state could be potentially probed to identity the 
molecule. 

In general, the spatial-light -modulator system of FIG. 4 
could be used to produce more general pulse shapes than a 
series of beats at a single resonance frequency. By using a 
more complicated pulse shape, several bonds present in a 
molecule can be coherently excited, and energy transferred 
between them in a coherent fashion. Because each molecule 
has some difference between the bonds presents and their 
relative orientation (and therefore the coupling in the quan- 
tum wave functions between them), pulses can be shaped that 
will selectively transfer energy between the states for a par- 
ticular molecule, and not be selective for other similar mol- 
ecules. In this way, the emission of stimulated Raman scat- 
tering or another coherent scattering process can be made 
more specific than just every molecule possessing a bond of a 
particular energy. With an automatically controlled pulse 
shaper, such as those based on spatial light modulators, feed- 
back can be employed where the computer can test various 
pulse shapes, measure the resulting emitted light temporal 
signal using the demodulator, and progressive reshape the 
pulse to optimize the signal from the target molecule and 
exclude other molecules. Once a useful temporal field shape 
for stimulating a molecule has been found, it can be stored in 
a database and later used for identifying that molecule in the 
future. 

2. Reference Generator. 

The reference generator takes a portion of the signal pro- 
duced by the oscillator and converts it to a reference signal. 
This reference signal is later used to demodulate the sample 
signal. The reference generator is a nonlinear process that 
converts light in the illumination bandwidth to light in the 
sample’s emission bandwidth, so that interference can occur 
between them. This nonlinear process may or may not be 
CARS or CSRS. 

A common implementation of the Reference Generator 
would be to focus the oscillator excitation into a sample of the 
same molecular species that one wishes to image. The refer- 
ence signal should then be very similar to the same molecular 
species contained in the sample. This is because they are the 
same molecule, illuminated by nearly identical pulses, con- 
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verting them to the output signal using resonant CARS or 
CSRS. The magnitude of the cross-correlation between these 
two signals should be great because of their similarity. In 
addition, if there is variability of oscillator output, having the 
5 reference generator and sample contain the same substance 
will respond in similar ways, so that the cross-correlation 
signal can remain high despite fluctuations in the oscillator. 
The benefit of using the same molecule in the reference gen- 
erator is that it is the molecule’s signal that acts as its own 
to “fingerprint” with which the cross-correlation can use to rec- 
ognize the molecule in the sample. If more selective excita- 
tion processes than CARS are used, then using the same 
molecule in both reference and sample will help ensure that a 
reference signal can be generated for a given excitation field. 
15 A nonresonant nonlinearity can be used as the reference 
generator as long as the peak power of the oscillator signal 
can excite a sufficient quantity of reference signal. High peak 
power can be maintained by not chirping the oscillator signal 
that is sent to the reference generator, while sending a rela- 
20 tively low peak power signal to the sample. Nonresonant 
CARS or CSRS can be used to generate an anti-Stokes or 
Stokes signal, respectively, from a medium that does not 
necessarily have a resonance at the frequency of the target 
molecule. The benefit of this is that the medium may not have 
25 to be changed each time a different molecular species is 
scanned for, because otherw ise a medium with a resonance at 
that wavelength would need to be chosen. Also, this species 
can act as a standard signal source against which the return 
signals from many samples can be compared. The nonreso- 
30 nant CARS can be implemented by focusing the excitation 
radiation into a sample of liquid that produces a CARS/CSRS 
signal in the same frequency band as that generated from the 
sample. For example, benzene will generate a CARS anti- 
Stokes signal in the 3000-3100 cm -1 frequency band. 

35 Continuum generation is another type of nonresonant non- 

linear process that can be used in the reference generator. A 
sufficiently high peak power pulse is focused into a medium, 
where it excites a broad bandwidth of frequencies to be pro- 
duced. If the produced frequencies overlap the emission fre- 
40 quency band produced in the sample, this portion of the 
continuum can act as a reference signal. The generated con- 
tinuum must be created by a mechanism that is sufficiently 
stable to not be overly sensitive to fluctuations in oscillator 
intensity. An unstable reference signal will result in noise in 
45 the cross -correlation signal. The benefit of continuum gen- 
eration is that is likely to create a broad bandwidth that will 
have signal in the emission bandwidth of the sample, so that 
the continuum need only be filtered for the needed frequency 
band. Also, if the oscillator employs continuum generation, it 
50 may already generate light within the emission bandwidth 
that can be used as a reference generator, eliminating a sepa- 
rate nonlinear process in the reference generator step. Some 
examples of materials used for continuum generation mate- 
rials are optical glass, fused silica, calcium fluoride, sapphire, 
55 ethylene glycol, water, high numerical aperture optical fibers, 
photonic crystal optical fibers, microstructured optical fibers, 
dispersion-shifted optical fibers, and gas cells (e.g. cells filled 
with helium, argon, or nitrogen). 

Other candidate processes for nonresonant nonlinear ref- 
60 erence generation include second and higher harmonic gen- 
eration, stimulated Raman scattering, sum and difference 
frequency generation, optical parametric amplification, four- 
wave mixing, and self phase modulation. 

FIG. 5 shows an implementation of a reference signal 
65 generator 500. The concentration optics 503 are typically 
implemented as some combination of lenses and mirrors. The 
concentration optics may also require some combination of 
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frequency dispersive elements such as prisms, diffraction 
gratings, pulse shapers, and dispersive materials to prepare 
the temporal shape of the signal for nonlinear generation. 
Concentration in space and time may be necessary because 
the nonlinear processes are power sensitive, and depends on 
the strength of the nonlinear process. When the light 501 
enters the nonlinear medium 505, it undergoes conversion to 
a frequency band coinciding with the frequency band of the 
response signal from the sample. This nonlinear medium may 
be one of the media mentioned above, either a sample of a 
target molecule, a solvent, or a continuum generation 
medium. After exiting the nonlinear medium, the reference 
signal 509 is collected by reference collection relay optics 
507, where it is sent to the demodulator where it is combined 
with the sample signal. The collection relay optics are usually 
implemented as some combination of lenses and mirrors that 
collimate the reference field radiation. This reference field 
should be characterized to find its temporal structure by 
instruments such as Frequency Resolved Optical Gating, 
cross-correlation with another known signal, nonlinear sono- 
grams, or nonlinear autocorrelations/cross -correlations. 

3. Microscope. 

The microscope delivers the excitation radiation from the 
oscillator to the sample and collects the resulting sample 
emission. Inside the sample, a coherent nonlinear process 
such as CARS or CSRS takes place that, in the presence of a 
molecule of interest, will emit the sample field in response to 
the excitation field. The sample field is collected by the micro- 
scope and then propagated to the demodulator, where the 
sample field can be estimated from the measured cross-cor- 
relation between the sample field and reference field. 

Microscope systems can be differentiated by various 
implementation choices. They can either illuminate one (se- 
rial scanning) or many points (full field imaging) at a time on 
the tissue. If the pump and Stokes (anti-Stokes) beams of the 
excitation field are separated in frequency, they can be sent in 
either separate (non-collinear) or identical (collinear) angles 
into the sample. The temporal delay of the response radiation 
relative to the reference may or may not be used to range 
molecular constituents in the tissue. The response radiation 
can be collected in the forward scattering (forward CARS/ 
CSRS) or backward scattering (epi-CARS/CSRS) directions. 

The microscope measures spatially resolved molecular 
density by illuminating various points on the tissue with the 
oscillator field, collecting the emitted sample field, and 
recombining with the reference field in the demodulator. One 
point in the tissue may be illuminated at a given time, result- 
ing in serial or raster scanning of the molecular density 
through the tissue. Alternatively, a line or a complete plane of 
points may be illuminated, so that data may be acquired from 
many points in parallel. Illuminating and measuring the radia- 
tion from an entire plane of points is called full field imaging. 
At the time of this writing, lull field imaging is seldom used 
because it requires an array intensity detector such as a 
charge-coupled-device (CCD) to simultaneously measure the 
demodulated signals of all of the illuminated points. Unfor- 
tunately, as of this writing CCD arrays produce thermal dark 
noise at each pixel, and also have a relatively limited dynamic 
range of measurable intensity values. Demodulated interfer- 
ence signals often require very high dynamic range detection. 
It is conceivable that future CCD or other types of focal plane 
arrays (e.g. CMOS arrays) may overcome these limitations. 
Full field imaging also requires that the tissue be illuminated 
by larger amounts of power because measurable signal must 
be produced for an entire area rather than just one point. Since 
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this is more likely to result in tissue damage, full field imaging 
will probably be used when speed of acquisition is para- 
mount. 

FIG. 6 shows three examples of full field CARS configu- 
5 rations. The “beam delivery optics” are usually implemented 
as some combination of mirrors and lenses that deliver beams 
that illuminate a wide area or line on the sample. For all of 
these microscope configurations, beam delivery and collec- 
tion optics will typically utilize a microscope objective. The 
to “response field collection optics” 601 are similarly imple- 
mented as a combination of lenses and mirrors that relay the 
response field to the demodulator so that it may be recom- 
bined with the reference and detected. The noncollinear full 
field CARS 602 delivers the pump 603 and the Stokes 605 
1 5 beams (assuming they are separate) at separate angles, so that 
the response is separated in angle from the illumination by an 
angle given by phase-matching conditions. The collinear 
geometry 607 sends the pump and Stokes radiation in the 
same directions (or in a single field if they can not be sepa- 
20 rated) and collects the radiation in the same direction, which 
can be discriminated with a dichroic beamsplitter. The epi- 
CARS geometry 611 collects the backscattered radiation, 
usually through the same objective optics that the sample is 
illuminated through. The epi-CARS can be discriminated 
25 from the illumination with a dichroic beamsplitter or an inter- 
ference filter. 

If collinear CARS is used, where the pump and Stokes 
(anti-Stokes) beams from the excitation beams overlap, then 
the response field 609 will overlap the excitation beams. Then 
30 the response field frequency band should not completely 
overlap the excitation frequency band, so that a spectral filter 
may distinguish between the excitation and response fields. If 
non-collinear CARS is used, then the response beam can be 
sufficiently angularly separated (as determined by the phase- 
35 matching criterion) from the excitation radiation to be filtered 
by a spatial filter. However, in the case of non-phase-matched 
CARS, such as epi-CARS, the interaction CARS volume 
must be small enough so that the response is radiated effec- 
tively isotropically, so that spatial filtering is unavailable and 
40 spectral filtering should be used. 

The other, more commonly used alternative is serial point 
scanning. Serial scanning tightly focuses the oscillation sig- 
nal into the tissue to create a very small volume where peak 
power is maximized. The focusing is usually achieved using 
45 a microscope objective. This focus is then scanned through a 
3-D set of points in the tissue, and the sample signal gathered 
from each point is demodulated to produce a 3-D NIVI image. 
Since nonlinear processes are power sensitive, efficient 
CARS/CSRS occurs only at the focus. If the sample is small 
50 enough, the focus may be scanned through the sample by 
translating the sample in all three dimensions. However, it is 
not feasible to move large samples such as human subjects 
this way. In this case, the beam focus can be moved in the 
transverse direction by steering the beam, perhaps using gal- 
55 vanometer rotated mirrors, acousto-optic modulators, or 
translating the lens assembly. The depth may be scanned by 
mechanically adjusting the distance between the lens and the 
tissue, perhaps using a lead-screw translator and/or a piezo- 
electric transducer. Since the signal can be excited at only one 
60 point at a time, one can be sure that the resulting measured 
sample signal is due to the interference of emissions of radia- 
tion produced in that volume only. This can be a benefit in 
NIVI when high phase resolution is required because one is 
assured that any measured phase shifts are not due to inter- 
65 ference between molecules at disparate spatial locations . This 
higher phase resolution may be used to better differentiate 
between similar molecular species. Serial point scanning 
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typically utilizes a single photodetector or a small number of 
photodetectors at the demodulator, which has the benefit that 
the dark current of a single photodetector is usually less than 
that of an entire CCD array, and a single photodetector can 
also typically handle a higher dynamic range of measure- 
ments. 

FIG. 7 shows the geometry of translated serial-port scan- 
ning configurations. In all of these configurations, the sample 
is on a translator that moves the sample through the focus to 
form an image of the molecular density at various points. The 
translator could be a three-axis linear screw drive translator, 
or piezoelectric translator, or a combination of these. The 
beam delivery optics 701 focuses the pump 703 and Stokes 
(anti-Stokes) 705 beams onto the tissue at a point of interest, 
and the response field collection optics 707 gathers the gen- 
erated anti-Stokes (Stokes) field 709 from the tissue. In the 
non-collinear geometry 700, the pump and Stokes illuminate 
the point of interest at different angles, so that the anti-Stokes 
will emerge at a third angle given by phase-matching, so that 
the anti- Stokes is spatially separated from the illumination. In 
the collinear geometry 702, the illumination and response 
fields emerge overlapped, so that they must be discriminated 
by frequency (e.g. using a dichroic beamsplitter). Finally, in 
the epi-CARS geometry 704, the backscattering response 
fields are collected, often by the same optics through which 
the illumination is projected, and can be separated with a 
dichroic beamsplitter. 

Alternatively, the focus can be moved and the sample can 
be left stationary. This can be accomplished by a combination 
of tilting the illumination beams before it enters the beam 
delivery optics and/or translating the beam delivery optics 
around (FIG. 8). Translating the collection optics and/or tilt- 
ing the exiting response beam with a beam steering optics 801 
can capture the exiting response field. The beam steering 
optics changes the direction of the incoming beam 803. This 
can be implemented, for example, by a galvanometer- 
scanned rotating mirror. By changing the direction of the 
beam before it enters the beam delivery optics 809, the posi- 
tion of the focus in the sample can be changed. The translation 
for the beam delivery and collection optics 807 can be imple- 
mented with a piezoelectric and/or a linear screw-drive trans- 
lation stage. Translating the beam delivery and/or collection 
optics moves the focus with the optics through the sample. 
These two mechanisms can be combined to enable the three- 
dimensional translation of the focus through the sample. This 
configuration is especially convenient for the epi-CARS con- 
figuration, because the same beam steering optics and beam 
delivery optics 809 can be used to collect the response field 
805. For these configurations, a compensating delay may 
need to be incorporated into the demodulator because steer- 
ing the beam and/or translating the delivery and collection 
optics change the travel time for the illumination and 
response signal through the microscope section. 

All of the previous scanning modes, full field imaging, 
translated serial -point scanning, and beam-steering serial- 
port scanning, use the spatial location of the illumination 
beam to differentiate between response signal gathered from 
various points in the sample. This confinement method is the 
same as that used by multiphoton microscopy or CARS 
microscopy. Other technologies such as Optical Coherence 
Tomography and Optical Coherence Microscopy use tempo- 
ral ranging in addition to spatial confinement to further isolate 
the contributions of signal from various points in the sample. 
Because of the heterodyne nature of NIVI, this scanning 
mode is also available. It may be attractive for in vivo imaging 
because it will enable scanning tissue without translating the 
microscope objective or sample, and may scan faster because 


16 

there are mechanisms for scanning the temporal delay much 
faster than translating objective optics. The phase measure- 
ment capability of NIVI is useful for both temporal ranging 
and molecular species identification this way. 

5 FIG. 9 shows two configurations of NIVI 900 and 902 
utilizing temporal ranging. Temporal ranging is achieved by 
measuring the interference of the reference and response 
signals for various relative temporal delays. The temporal 
gating configuration superficially resembles the other epi- 
10 CARS configurations. However, unlike previous scanning 
modes, the depth-of-field of the focusing of the illumination 
from the beam delivery optics in the tissue is set to be long, 
because the temporal gating will be used to discriminate 
between molecular constituents within the depth-of-field. 
15 This temporal ranging utilizes epi-CARS/epi-CSRS because 
the backscattered response signal is collected from the 
sample by the response field connector so that the signal delay 
into the tissue can be timed. Because CARS/CSRS is not 
phase-matched for the backwards direction, the generation of 
20 a backscattered signal will only occur efficiently for small 
particles or edges of particles with a compatible resonance. 
The chief difference between this configuration and the other 
epi-CARS configurations is that the interference signal is 
scanned as a function of relative delay, and an interference 
25 signal maximum indicates the presence of a molecular spe- 
cies at a particular depth in the medium (corresponding to that 
time delay). In the full field configuration 900, an entire plane 
of points is interfered with the reference signal to produce 
simultaneous measurements of the cross-correlation signal 
30 for the entire plane 909 . By scanning the delay, the molecular 
density of various planes can be measured. For the beam- 
steered temporal gating setup 902, the beam is scanned 
through the sample by tilting the beam through the beam 
steering optics 903. Beam steering provides lateral displace- 
35 ment of the beam, and temporal gating provides depth infor- 
mation so that a three-dimensional volume is scanned. The 
sample can be translated laterally also to scan the beam. The 
implementation differences between the temporal gating con- 
figurations are in the choice of depth-of-field of the illumina- 
40 tion and collection optics, and in the demodulator. The 
demodulator must be designed to scan a sufficient temporal 
interval to capture the interference signal between the refer- 
ence and response signals. This temporal interval is typically 
between 500 microns to 10 mm. 

45 4. Demodulator. 

The purpose of the demodulator section is to decode the 
response signal from the sample so that it may be measured by 
relatively slow electronic equipment (slow compared to the 
oscillations of the electric field of the response signal). The 
50 magnitude of this demodulation signal will be related to the 
density of a molecular species of interest in the tissue. With 
knowledge of the molecular density at each point, a molecular 
density map, or NIVI image, can be presented to the user. This 
demodulation is implemented as the cross-correlation of the 
55 response signal relative to a generated reference signal. Out- 
lined below are various optical configurations that produce 
this cross-correlation signal. 

There are various design choices that are made when 
choosing a demodulator. First, one needs to know whether or 
60 not full field microscopy is used. Also, one must decide 
whether the cross-correlation will be measured one sample at 
a time, or many samples in a single instant. 

FIG. 10 shows an example of a cross -correlator that can be 
used to demodulate the full field CARS signal, measured for 
65 example as in FIG. 6. The response field is collected and 
relayed by the “response field collection optics” of FIG. 6 to 
the Response Field 1001 of FIG. 10. The reference field 1003 
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is produced by the Reference Generator and relayed to the 
cross -correlator. The reference is delayed relative to the 
response field, and they are mutually overlapped using a 
beamsplitter 1005. The “response field spectral filter” 1007 
filters the recombined field for only the frequency band that 
contains the response field bandwidth, which removes any 
remaining oscillator signal. The filter can be implemented by 
a combination of interference filters, color glass filters, dich- 
roic beamsplitters, or other frequency selective elements. The 
combined fields are imaged onto a focal plane array 1011 , 
such as a CCD, so that the CCD is the conjugate image plane 
of the sample plane. The imaging is achieved with the “relay 
imaging optics” 1009 which are some combination of lenses 
and mirrors. The intensity detected on the CCD corresponds 
to a single cross -correlation signal measurement collected 
from various points in the sample. The “variable attenuator” 
1013 is adjusted to maximize the use of the dynamic range of 
intensity measurements of the CCD. As the focus in the 
sample changes (e.g. by translating the sample), and/or the 
adjustable delay 1015 is changed, the entire cross -correlation 
signal can be measured, forming a three- -dimensional data set. 
The typical scan range length for high numerical aperture full 
field imaging will be up to 100 microns. The delay can be 
implemented as a mirror translated by a linear-screw-drive 
translation stage, or by piezoelectric actuation. This same 
cross-correlation can also be used for temporal ranging full 
field NIVI, as shown in FIG. 9, with the only difference being 
that the delay mechanism must be designed to scan a suffi- 
ciently long range of interest in the sample, typically from 
500 microns to 10 mm. 

The measurement of a cross-correlation signal for a serial- 
point scanning microscope is basically the same system as for 
full field imaging, except that a single photodetector can be 
used rather than an array photodetector. The adjustable delay, 
response field spectral filter, and relay imaging optics can all 
be implemented in similar ways to the full field case. The 
relay imaging optics will need only focus the combined 
response/reference signal onto the photodetector. In some 
cases, when only the magnitude of the cross-correlation sig- 
nal at its peak needs to be measured, it may be desirable to 
dither the adjustable delay with a piezoelectric transducer a 
fraction of a wavelength, so that the magnitude of the cross- 
correlation signal peak can be demodulation with a multiply- 
ing mixer and low-pass filter (an electronic heterodyne 
demodulator). This configuration will be covered in more 
detail in part five. To measure the entire cross -correlation, the 
adjustable delay will be scanned over various time delays and 
the photodetector intensity signal measured. 

In FIGS. 10 and 11, a variable attenuator 1013 is used to 
adjust the intensity of the signal that reaches the photodetec- 
tor 1101 so that its dynamic range is not exceeded. Variable 
attenuators can be implemented with liquid crystal shutters, 
neutral density filter wheels, or rotating polarizers. The spa- 
tial filter 1103 is used on the response signal to filter out 
spatial inhomogeneities in the response signal that could 
reduce the depth of modulation at the photodetector. A spatial 
filter would typically consist of a telescope of two converging 
lenses, with a pinhole in the focal plane between the lenses to 
filter out power around the main diffraction focus. 

When the oscillator produces pulses of a low repetition 
rate, so that for example high peak power can be employed, 
one may want to measure multiple samples of the cross - 
correlation signal simultaneously. This can be achieved by the 
configurations in FIG. 12. These configurations have the 
advantage that, with a sufficient amount of pulse power and 
number of photodetector array samples, the cross -correlation 
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signal can be measured using a single pulse. Low repetition 
rate illumination can keep the peak power high while average 
power remains low. 

The “linear photodetector array cross-correlator” 1201 
5 expands the response and reference signals, and interferes 
them with an angle between the two beams. The beam expan- 
sion is achieved by, for example, a pair of converging achro- 
matic lenses 1203 arranged in a telescope configuration. The 
points on the wavefront where the two signals combine will 
to have various time delays between them. A cylindrical lens 
1205 then focuses the beams into a line image on a linear 
photodiode array 1207. This concentrates the signal to adjust 
for the narrow height of the detector pixels. Each intensity 
sample on the linear photodetector array indicates a sample of 
15 cross -correlation of the two fields with various relative time 
delays, with a constant intensity signal added. The recorder 
would then read the linear CCD signal so that a complete 
NIVI image can be assembled. 

The primary disadvantage of this system is that the modu- 
20 lation of the intensity signal on the linear photodetector array 
is rather low, and the dynamic range of typical CCDs is 
likewise low, so that the signal cannot be measured to a high 
signal-to -noise ratio. To combat this, one can utilize the Fou- 
rier- transform cross-correlator 1209. Rather than directly 
25 measuring samples of the cross-correlation, it interferes the 
two signals and measures their spectral decomposition. The 
beamsplitter 1005 combines the response and reference sig- 
nals, which are filtered for the response bandwidth by utiliz- 
ing a response field spectral filter 1007. This signal is then 
30 filtered by a frequency dispersive element 1211 such as a 
diffraction grating that scatters each frequency to a different 
angle. The power of each frequency is then focused to a 
different pixel on the linear photodetector array 1215 by using 
a focusing element 1213 (typically a combination of lenses or 
35 mirrors). The samples of the intensity measured on the pho- 
todetector indicate the real part of the Fourier transform of the 
cross-correlation signal. The recorder may compute an 
inverse discrete Fourier transform of the linear photodetector 
array intensity samples to recover the cross-correlation. The 
40 Fourier transform already performs the necessary Hilbert 
transform to infer the imaginary part of the cross -correlation 
signal. However, because it is an intensity measurement, the 
cross-correlation signal, a time-reversed version, and the 
autocorrelation of the response signal are superimposed in the 
45 intensity signal . By choosing the adjustable delay longer than 
the length of the cross -correlation temporal signal, the recon- 
struction of these signal components will not overlap in the 
time domain, and so the recorder device may distinguish the 
cross -correlation from its mirror image and the autocorrela- 
50 tion. 

The adjustable delay may also be dithered a small amount 
to introduce a phase shift into each measured frequency com- 
ponent at the linear photodetector array. By utilizing at least 
three linearly independent phase shifts, the phase of each 
55 Fourier component can be known and the time domain cro ss - 
correlation computed using an inverse discrete Fourier trans- 
form. However, then the ability to measure the cross-correla- 
tion of a single pulse will be lost. Another possibility is to use 
a linear photodetector array with several rows (at least three) 
60 rather than just one row of pixels. By tilting one of the wave- 
fronts slightly vertically with respect to the other, a small 
phase shift can be introduced in the measured cross-correla- 
tion signal between rows on the linear photodetector. With a 
sufficiently large phase shift, the complex amplitude of each 
65 Fourier component can be inferred from the intensity samples 
from each column on the linear photodetector by utilizing the 
discrete Fourier transform. Since all of the time delayed 
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cross -correlation samples are measured simultaneously, this 
is a single-pulse measurement. 

To minimize the effect of dark charge built up in the linear 
CCD detector, it is best to discard the charge as soon as 
possible before the pulse arrives, and read the CCD as soon as 
possible after the pulse arrives. The speed of readout should 
be as fast as possible given limitations in accuracy due to the 
readout noise. 

5. Record r. 

The recorder accumulates the samples of the cross-corre- 
lation signal gathered from various points in the sample, 
processes this data, and presents these as a human-interpret- 
able image. It is usually implemented as a data acquisition 
digital computer with some method of automatic control of 
the adjustable mechanisms (such as delay lines, piezoelectric 
actuators, and galvanometer mirrors), analog-to -digital con- 
version, and some form of image output device such as a 
screen or printer. It may also have control of the oscillator 
itself, to automatically tune the wavelengths or bandwidths of 
output or control the rate or timing of pulse output. It may also 
control the delay lines or spatial light modulators in the pulse 
shaping mechanisms of FIG. 4. 

Typically, the recorder will scan the illumination through 
the sample and/or the adjustable time delay and measure the 
cross-correlation signal. Unless the Fourier-Transform cross- 
correlator is used, the intensity samples represent the cross- 
correlation signal with a constant level added that could be 
discarded. The magnitude of the cross -correlation signal indi- 
cates the presence of the target molecular species. For CARS/ 
CSRS the magnitude of the cross -correlation signal is related 
to the second power of the molecular density of that species. 
By demodulating the magnitude of the cross -correlation sig- 
nal, the molecular density can be estimated for the points on 
the cross -correlation signal for the areas from which the 
response signal was collected to form that cross -correlation 
signal. Both spatial confinement and temporal ranging can be 
used to differentiate between the signals due to molecular 
densities at different locations in the sample. 

FIG. 13 documents the relationship between the elements 
of the recorder 1300. For serial -scanning or temporal-scan- 
ning configurations, the instructions from the operator are 
entered via a human interface device 1304 into the digital 
computer 1301. The digital computer controls the scanned 
beam position in the sample using the galvanometer mirror 
angles of the galvanometer scanned mirrors 1303, the oscil- 
lator 1302, the adjustable delay line 1309, the position of the 
microscope objective position 13 05, and/or the position of the 
sample translator 1307 to scan the illumination through the 
sample. In general only a subset of these need be controlled to 
scan the three-dimensional volume of the tissue. If full field 
imaging is used, usually only the depth and/or the delay line 
need be scanned. The cross-correlation signal as measured by 
the intensity is converted to a voltage by the photodetector 
1311, which is in turn converted to a digital sample by the 
analog-to-digital converter 1313. The digital computer varies 
the delay line and/or reads various pixels from the photode- 
tector (if it is an array detector) to determine the cross-corre- 
lation signal. If a Fourier- transform cross-correlator is used, 
the computer will need to compute the inverse discrete Fou- 
rier transform of the signal to find its cross-correlation from 
the intensity samples. The computer then associates the mag- 
nitude of the cross-correlation signal with the molecular den- 
sity, and assembles these magnitudes into a density map of the 
sample. This density map is stored in the storage device 1312 
presented on the visual display 1315 and/or made into a 
physical representation with a printer or stereolithography 
device 1317. 
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To aid in the measurement of the magnitude of the cross - 
correlation signal, the configuration 1400 of FIG. 14 is sug- 
gested. 

Because the digital computer 1301 can record only rela- 
5 tively slow signals, a dither oscillator 1401 may be used to add 
a small high-frequency dither signal into the adjustable delay 
line 1309 that produces a periodic perturbation of the delay in 
the signal of a magnitude usually of less than one wavelength. 
This same signal is multiplied by the received photodetector 
to 1403 interference signal to demodulate it, and is low pass 
filtered through the low pass filters 1405 to remove harmonics 
of the dither frequency. Both the dither and its quadrature 
signal are demodulated, because these correspond to the real 
and imaginary parts of the complex cross-correlation ampli- 
15 tude. The analog-to-digital converter 1313 will then be 
directly measuring a quantity proportional to the complex 
magnitude of the cross-correlation for the delay line position. 
The computer can utilize the digitized real and imaginary 
cross -correlation components to display the amplitude and 
20 phase of the cross -correlation signal. The amplitude will cor- 
respond to the magnitude of the reflection, and the phase will 
correspond to the Doppler shift of the reflection. If a very fast 
dither signal is desired (above 10 kHz), an electro-optic 
modulator or acoustical-optic modulator can be placed in the 
25 adjustable delay line system to rapidly modulate the delay a 
small amount. 

The magnitude of the cross-correlation signal depends 
both on the density of molecules available to produce the 
response signal, and on the temporal structure of the signals 
30 radiated from the molecules. The temporal signal produced 
by the molecules can be inferred from the cross -correlation 
signal. If f(t) is the temporal response signal produced by a 
molecule, and g(t) is the known reference signal, the mea- 
sured cross-correlation T(t) is given by: 


r(r) = J f(t)g{t-r)dt 
40 

If F(oo), G(go), and T((o) are the Fourier transforms of f(t), 
g(t), and T(x) respectively, then r(to)=F((n)G(-0L)). The func- 
tion f(oo) can be computed from the Fourier transform of the 
cross -correlation and G(co) can be computed from the mea- 
45 sured reference signal. The Fourier transform F (oo) of tempo- 
ral signal f(t) can then be estimated by F(w)=T((b)G(-ol))7( 
G(-cd) 2 +N(go) 2 ) (change to equation), where N(oo) is an esti- 
mated power- spectral-density of the noise. An inverse Fourier 
50 transform of F(oo) yields f(t). 

The recorder can perform this computation and thus 
recover f(t) for various molecular species. An unknown mol- 
ecule may be identified by comparing its temporal signal to 
known molecular signals. In addition, by looking at the fre- 
55 quencies co of phase arg F(oo) where the phase changed rap- 
idly, the resonance frequencies of the molecules of interest 
can be determined to high precision. A library of the temporal 
responses and resonance frequencies of various molecules for 
CARS excitation may be built up and used to identify 
60 unknown molecules in vivo . We also note that the same cro ss - 
correlation measurement configurations can be used to mea- 
sure the OCT backscattered signal due to the scattering of the 
excitation radiation. To do this, the excitation and response 
bands should be separated with a dichroic beamsplitter, and 
65 the cross -correlation setups implemented separately to elimi- 
nate photon noise from the excitation band being measured in 
the response band. Different cross-correlation setups may be 
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used for each configuration, e.g. a single-pulse measurement 
system for the conventional excitation and a time delay cross- 
correlator as in FIG. 11 . 

Embodiment 

Randomness introduced into the measurement process 
may produce false positive indications of molecular densities, 
or obscure weak concentrations of molecular densities. This 
randomness has three sources: fluctuations in the oscillator, 
vibrations and air currents, and noise introduced by the pho- 
todetector. 

Oscillator fluctuations are the hardest to characterize 
because the feedback mechanism of laser sources can pro- 
duce large variations in the frequency, bandwidth, and output 
power of the pulses even for small perturbations of the oscil- 
lator. In addition, nonlinear processes such as self-phase- 
modulation, self-focusing (Kerr lensing), or continuum gen- 
eration may also affect the pulse in ways that are sensitive to 
the power in the pulse. Therefore, it is desirable to keep 
fluctuations in pulse energy below a few percent, and keep 
frequency and bandwidth fluctuations under one percent. 

Well designed state-of-the-art Ti-sapphire oscillators and 
regenerative amplifiers can meet these specifications if they 
are aligned and maintained within the guidelines set up by the 
manufacturer. It is important that the internals of the appara- 
tus be shielded, as well as the paths of the beams, as much as 
possible from air currents and vibrations. This can be 
achieved by placing the apparatus on a platform with air 
flotation or shock-mounted legs. The oscillator apparatus 
should be enclosed in a rigid case with baffles that minimize 
the opening surface area to that required for the beams to enter 
and leave the case. It is best to use laser-diode or laser-diode- 
pumped solid-state sources to pump the oscillator to mini- 
mize pump fluctuations, but a well controlled ion laser is also 
usable. 

A dispersion compensated mirror Ti-sapphire system can 
produce 200 nm or more bandwidth centered at 800 nm and is 
also suitable for the chirped-CARS or CARS utilizing a pulse 
shaper. 

Another possibility is to employ a femtosecond diode- 
pumped erbium or ytterbium fiber laser. Fiber-based con- 
tinuum generation sources of nanojoule energy have been 
demonstrated, and this source may be relatively inexpensive 
and compact. A pulse selector (“pulse picker”) utilizing a 
Pockels cell or a fiber-based electro-optic modulator or elec- 
tro-absorption modulator may be used to lower the average 
power while keeping the peak power high. A fiber amplifier 
can be used to increase the energy of the pulses. A fiber-based 
amplifier may provide a more portable, stable source. 

If continuum generation is used, a source employing self- 
phase-modulation is preferable because of the short interac- 
tion length of the nonlinearity. For example, a microjoule 
energy pulse will broaden when focused into a fused silica, 
sapphire, calcium fluoride, or quartz medium. Self-focusing 
can be used to increase the peak power and minimize the 
interaction length. Self-phase-modulation in an optical fiber 
may be sufficient if it occurs over a short length of fiber. The 
fiber should be single mode in the bandwidth utilized to 
minimize nonlinear coupling of multiple modes. Higher peak 
power allows a stronger nonlinearity to be used and therefore 
minimizes the interaction length. 

Photodetector noise is the other source of noise. There are 
two sources of photodetector noise: photon noise and thermal 
noise. Thermal noise may be eliminated by using a cooled 
detector, or by minimizing the integration time. If a line 
camera is used, the charge on the array can be expelled before 
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the pulse is received and the signal can be read out immedi- 
ately after the pulse is received. Also, the reference power 
should be balanced with the power received from the sample 
so that the dynamic range of the detector is utilized. Photon 

5 noise is fundamental to photodetection processes. Its effect 
can be minimized in several ways: using more illumination 
power to produce a larger CARS signal. Also, using the Fou- 
rier-transform cross -correlator minimizes the photon noise 
effect by increasing the modulation when single pulse illu- 

10 mination is used. 

Equipment Brand Names: 

Coherent Laser (of Coherent, Inc.) 10W Verdi pumps the 

Kapteyn-Mumane Laboratories, LLC, titanium: sapphire 
seed (pump) laser 

1 5 Which pump s the 

Coherent Regenerative Amplifier (RegA 9050) 

And uses the 

Coherent Stretcher/Compressor. 

20 A portion of the amplified light is sent to a 

Coherent Optical Parametric Amplifier (OPA9450) 

Which produces the second optical beam. 

While we refer to a microscope (Leica) for beam delivery 
and holding the sample, any standard optics could be used. 

25 We are currently using a computer controlled 3 -axis stage 
(Newport Corp.) to translate our samples under the imaging 
beams. 
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50 The invention claimed is: 

1. A method of examining a sample, comprising: 
exposing a non-linear reference medium to a first set of 
electromagnetic radiation generated in an interferomet- 
ric system to form a second set of electromagnetic radia- 
55 tion scattered from the non-linear reference medium; 

exposing a sample to a third set of electromagnetic radia- 
tion to form a fourth set of electromagnetic radiation 
scattered from the sample; and interfering the second set 
of electromagnetic radiation and the fourth set of elec- 
60 tromagnetic radiation; and 

detecting fourth photos in the fourth set of electromagnetic 
radiation; 

wherein the detecting comprises interfering; 
wherein the fourth photons are anti-Stokes photons or 
65 Stokes photons; 

wherein the first set and the third set of electromagnetic 
radiation are generated from a source; 
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at least a portion of the second set of electromagnetic 
radiation is of a frequency different from that of the first 
set of electromagnetic radiation; and 
at least a portion of the fourth set of electromagnetic radia- 
tion is of a frequency different from that of the third set 5 
of electromagnetic radiation. 

2. The method of claim 1, wherein the fourth photons are 
anti-Stokes photons. 

3. The method of claim 1, wherein the fourth photons are 

Stokes photons. 10 

4. A method of forming an image of a sample, comprising: 
exposing a non-linear reference medium to a first set of 

electromagnetic radiation generated in an interferomet- 
ric system to form a second set of electromagnetic radia- 
tion scattered from the non-linear reference medium; 15 
exposing a sample to a third set of electromagnetic radia- 
tion to form a fourth set of electromagnetic radiation 
scattered from the sample; 

forming a digital data set corresponding to the sample; and 
converting the data set into an image; 20 

wherein the forming of the data set comprises interfering 
the second set of electromagnetic radiation and the 
fourth set of electromagnetic radiation; 
the first set and the third set of electromagnetic radiation 
are generated from a source; 25 

at least a portion of the second set of electromagnetic 
radiation is of a frequency different from that of the first 
set of electromagnetic radiation; and 
at least a portion of the fourth set of electromagnetic radia- 
tion is of a frequency different from that of the third set 30 
of electromagnetic radiation. 

5. The method of claim 4, wherein the fourth photons are 
anti-Stokes photons. 

6. The method of claim 1, wherein the electromagnetic 
radiation is in the frequency range of infra-red to ultraviolet 35 
light. 

7. The method of claim 4, wherein the electromagnetic 
radiation is in the frequency range of infra-red to ultraviolet 
light. 

8. The method of claim 1, wherein the examining of the 40 
sample is by optical coherence tomography. 

9. The method of claim 4, wherein the forming of the data 
set is by optical coherence tomography. 

10. The method of claim 1, wherein the sample is selected 
from the group consisting of a tissue sample, a single cell, and 45 
a patient. 

11. The method of claim 4, wherein the sample is selected 
from the group consisting of a tissue sample, a single cell, and 
a patient. 

12 . The method of claim 1 0, wherein the sample is a human 50 
patient. 

13 . The method of claim 11, wherein the sample is a human 
patient. 

14. In a method of forming an image using an optical 

coherence tomography, apparatus, the method including 55 
exposing a sample or patient to electromagnetic radiation, 
collecting scattered electromagnetic radiation, and forming 
an image from the collected electromagnetic radiation includ- 
ing interfering the collected electromagnetic radiation with 
reference electromagnetic radiation, the improvement com- 60 
prising the reference electromagnetic radiation being scat- 
tered from a non-linear reference sample, and the wavelength 
of the collected electromagnetic radiation being different 
from that of the electromagnetic radiation that the sample or 
patient is exposed to. 65 

15. In an optical coherence tomography device, including 
an electromagnetic radiation source for generating reference 
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electromagnetic radiation and sample electromagnetic radia- 
tion, an optical delay line, a scanner, and a electromagnetic 
radiation detector, the improvement comprising a reference 
holder and optics adapted for exposing the reference electro- 
magnetic radiation to a non-linear reference medium before 
interfering the reference electromagnetic radiation with the 
sample electromagnetic radiation. 

16. A device for examining a sample, comprising: 
an oscillator, 

a reference generator, including non-linear reference 
medium, a optically coupled to the oscillator, 
a sample illuminator, optically coupled to the oscillator, 
an interferometric demodulator, optically coupled to the 
reference generator and the sample illuminator, 
a recorder, coupled to the demodulator, and 
a frequency- selecting element that ensures that the light 
that illuminates the sample is excluded from the modu- 
lator. 

17. The device of claim 16, wherein the oscillator com- 
prises a laser. 

18. The device of claim 17, wherein the oscillator further 
comprises a pulse shaper or a chirper. 

19. The device of claim 16, further comprising a scanner 
for scanning a sample, coupled to the sample illuminator. 

20. A method of examining a sample, comprising: 
exposing sample to a first set of electromagnetic radiation 

generated in an interferometric system to form a second 
set of electromagnetic radiation non-linearly scattered 
from the sample; and 

interfering the second set of electromagnetic radiation with 
a third set of electromagnetic radiation; 
wherein the third set of electromagnetic radiation is coherent 
with the first set of electromagnetic radiation; 

at least a first portion of the second set of electromagnetic 
radiation is of a frequency different from that of the first 
set of electromagnetic radiation; and 
at least a portion of the third set of electromagnetic radia- 
tion is of the same frequency as the first portion of the 
second set of electromagnetic radiation. 

21. The method of claim 20, further comprising detecting 
second photons in the second set of electromagnetic radia- 
tion; wherein the detecting comprises the interfering. 

22. The method of claim 21, wherein the second photons 
are anti- Stokes photons. 

23 . A method of forming an image of a sample, comprising: 
exposing a sample to a first set of electromagnetic radiation 

generated in an interferometric system to form a second 
set of electromagnetic radiation non-linearly scattered 
from the sample; 

forming a digital data set corresponding to the sample; and 
converting the data set into an image; 
wherein the forming of the image comprises interfering the 
second set of electromagnetic radiation and a third set of 
electromagnetic radiation; 

wherein the third set of electromagnetic radiation is phase- 
coherent with the first set of electromagnetic radiation; 
at least a first portion of the second set of electromagnetic 
radiation is of a frequency different from that of the first 
set of electromagnetic radiation; and 
at least a portion of the third set of electromagnetic radia- 
tion is of the same frequency as the first portion of the 
second set of electromagnetic radiation. 

24. The method of claim 23, further comprising detecting 
second photons in the second set of electromagnetic radia- 
tion; wherein the detecting comprises the interfering, and the 
forming of the image comprises the detecting. 
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25. The method of claim 23, wherein the second photons 
are anti- Stokes photons. 

26. The method of claim 23, wherein the second photons 
are Stokes photons. 

27. The method of claim 20, wherein the electromagnetic 
radiation is in the frequency range of infra-red to ultraviolet 
light. 

28. The method of claim 23, wherein the electromagnetic 
radiation is in the frequency range of infra-red to ultraviolet 
light. 

29. The method of claim 20, wherein the examining of the 
sample is by optical coherence tomography. 

30. The method of claim 23, wherein the forming of the 
image is by optical coherence tomography. 

31 . The method of claim 20, wherein the sample is selected 
from the group consisting of a tissue sample, a single cell, and 
a patient. 

32. The method of claim 23, wherein the sample is selected 
from the group consisting of a tissue sample, a single cell, and 
a patient. 

33 . The method of claim 31, wherein the sample is a human 
patient. 
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34. The method of claim 32, wherein the sample is a human 
patient. 

35. A method of examining a sample, comprising: 
exposing a non-linear reference medium to a first set of 

5 electromagnetic radiation generated in an interferomet- 
ric system to form a second set of electromagnetic radia- 
tion inelastically scattered from the non-linear refer- 
ence; 

exposing a sample to a third set of electromagnetic radia- 
10 tion to form a fourth set of electromagnetic radiation 
scattered from the sample; and 
interfering the second set of electromagnetic radiation and 
the fourth set of electromagnetic radiation; 
wherein the first set and the third set of electromagnetic 
15 radiation are generated from a source; 

at least a portion of the second set of electromagnetic 
radiation is of a frequency different from that of the first 
set of electromagnetic radiation; and 
at least a portion of the fourth set of electromagnetic radia- 
20 tion is of a frequency different from that of the third set 
of electromagnetic radiation. 
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